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Wigner Measures and Quantum Control 


AHHOTail,HH 

We study several examples from quantum control theory (see mm in the framework 
of Wigner functions and measures for infinite dimensional open quantum systems. An 
axiomatic definition of coherent quantum feedback is proposed within this setting. 


The representation of the states of quantum systems in terms of Wigner 
measures allows for useful analogies with classical system dynamics; it is similar 
to the representation of the states of classical Hamiltonian systems in terms of 
probability measures on the phase space. The passage to a description of the state 
of a subsystem of some larger quantum system is in both cases implemented by 
means of projective operations, because the phase space of the classical analogue of 
the ambient quantum system, being the union of some subsystems, is the Cartesian 
product of the phase spaces of the classical analogues of these subsystems. In the 
more familar setting of a finite dimensional phase space we can consider, instead of 
Wigner measures, their densities with respect to the Liouville measure, which are 
classical Wigner functions. However, on an infinite dimensional phase space, there 
exists no Liouville measure, i.e., a Borel cr-additive a-ffihte locally finite measure 
invariant with respect to symplectic transformations (this is a consequence of a 
well known theorem of Weil). In this case, we can either directly apply the Wigner 
measure, or introduce some “sufficiently nice” measure instead of the Lebesgue 
measure. In the case of a linear phase space, we can use a Gaussian measure as the 
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“sufficiently nice” measure, as is done in the so-called white noise analysis. After 
this, it becomes possible again to replace Wigner measures by “Wigner functions”, 
i.e., by their densities with respect to the new measure. We shall consider Wigner 
measures and their densities in parallel. 

Section 1, which is of independent interest, considers the properties of Wigner 
measures and functions; some of the results of this section can be regarded as an 
extension of results of [4] to Wigner measures. In section 2, an equation describing 
the evolution of the Wigner functions of quantum systems obtained by quantizing 
Hamiltonian systems with infinite dimensional phase space is given; this equation 
is obtained as a consequence of a similar equation for the evolution of a Wigner 
measure (see [5]). (A Wigner measure is a signed cylindrical measure, and it would 
be interesting to estimate its variation and find countable additivity conditions; 
however, we do not discuss these issues here.) 

The final section considers the evolution of the Wigner measures and functions 
of subsystems of quantum systems. Here, models of control of quantum systems 
are discussed and an axiomatic definition of coherent quantum feedback is given, 
which, as far as we know, has not been explicitly introduced in the literature 
so far. We consider largely algebraic aspects of the theory, omitting analytical 
assumptions. 

1. WIGNER MEASURES AND FUNCTIONS. 

This section discusses properties of Wigner measures and their densities with 
respect to fixed measures on a classical phase space, that is, Wigner functions 
(precise definitions are given below). Let E := Q x P be the phase space of 
a Hamiltonian system, where Q and P are real locally convex spaces (LCSs), 
P = Q*, and Q = P* (given an LCS X , we denote by X * its dual endowed with 
a locally convex topology consistent with the duality between X and X * ); then 
E* = P xQ. Suppose also that (•,•): PxQ i—> M is the bilinear form of the duality 
between P and Q. Then the linear mapping J : E 3 (q : p) i—» (p, q) G E * is an 
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isomorphism, and we identify h £ E with Jh £ E*. In particular, for each h £ E, 
the symbol h denotes the pseudodifferential operator on jC 2 (Q,p) whose Weyl 
symbol] is the function Jh £ E*. By p we denote the P-cylindrical (Gaussian) 
measure on j£ 2 (Q : p) whose Fourier transform : P —>■ M. is determined by 
&n(p) := exp (—| (p, P^p)), where : P —> Q is a continuous linear mapping 
such that (p, Bpp) > 0, for p ^ 0. By v we denote a Q-cylindrical measure on P 
whose Fourier transform dp : Q —» R is defined by $ v (q) := exp (—| (P*g,g)). 
In what follows, we assume that all LCSs are Hilbert, although the main results 
can be extended to the general case. We identify the space Q with Q* and P 
with P*, so that B* = B fl and > 0; note also that p and v are cr-additive 
if the operator B M is nuclear. The Weyl operator W(h) generated by an element 
h £ Eis defined by W(h) := e~ lh . The Weyl function corresponding to a density 
operator T is the function VFr(h) : E R defined by Wr{h) := tr (TW(h)) (see 
[4]); it does not depend on p. 

Definition 1 ([5]). The Wigner measure corresponding to a density operator 
T is an FT-cylindrical measure Wt on E determined by the relation 



e i (^)+^i))w T (dq 1 ,dp 1 ) 

p 


= y^T(h)(q 2 ,p 2 )- 


In other words, Wt is the (inverse) Fourier transform of the function Wr(h). 
Therefore we have W T (dq,dp ) = f Q fpWT(h)(q 2 ,p 2 )F ExE {dq 2 ,dp 2 ,dq,dp), 

where F ExE is the Hamiltonian Feynman psuedomeasure E x E. 

The Feynman pseudomeasure F E hou a Hilbert space is a distribution (in the 
sense of the theory of Sobolev-Schwartz generalized functions) on /C, on a Hilbert 
space is a distribution (in the sense of the theory of Sobolev-Schwartz generalized 
functions) on JC. It is convenient to specify such a functional F E . as well as an 
ordinary measure, in terms of its Fourier transform Fjc : /C 3 z i—>■ Fjc((f> z ) £ C, 
where p z : /C —>■ C is defined by ip z (x) := e l ^ z ' x \ 

4 The definition of a pseudodifferential operator F on C 2 (Q,J) with symbol F can be found in [5]. 
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HKL = E = QxP and Fjc(q,p) = e 1 ^, then F is said to be a 
Hamiltonian Feynman pseudomeasure; it is convenient for defining the Fourier 
transform that on functions given on infinite dimensional spaces and maps them to 
measures. Actually, the Hilbert space structure is not important here; a Feynman 
pseudomeasure, as well as a Gaussian measure, can be defined on any LCS; 
in particular, a Hamiltonian Feynman pseudomeasure can be defined on any 
symplectic LCS (additional information is contained in [3, 9, 11 ]). 

Proposition 1 (see [5]). If G is the Weyl symbol of a pseudodifferential 
operator on /^(Q, p), then f p f ( j G(q,p)Wr(dq, dp) = tr ^ TG^j . 

This proposition can also be used as a definition (cf. [4, Definition 3], where it 
is, however, assumed that dim Q = dim P < oc and, for this reason, only Wigner 
function, rather than measures, are considered). 

Definition 2. The density of the Wigner measure Wt with respect to q on 
Q x P (if this density exists) is called the Wigner 77 -function (if dim Q = dim P < 
00 and 77 is a Lebesgue measure on Q x P, to then the Wigner 77 -function is the 
classical Wigner function). 

In what follows, we assume that 77 = p <g> is, but refer to the Wigner p 0 v- 
function simply as the Wigner function. 

Corollary 1. If the assumptions of Proposition 1 hold, then 



/ / G(q,p)$ T {q,p)n® v(dq,dp) = tr(TG). 

Jp Jq 

Proposition 2. The following relation holds: 



K (^2 T/P2) + (<72 , W 1 92 



xe 2 



the Gaussian measure p <g) v (see [S] and the references therein). The relations 
given above and those similar to them can be obtained by using the following 
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heuristic rule. First, we write the corresponding formulas for the case where 
dimQ < oo, replacing Gaussian measures by their densities with respect to 
Lebesgue (=Liouville) measures on the spaces Q n Q x P; in turn, these formulas 
are obtained by using the standard isomorphisms between the spaces of functions 
square integrable with respect to the Lebesgue measure and the spaces of functions 
square integrable with respect to the Gaussian measures. After this, we pass to 
the infinite dimensional case, for which purpose we replace the Gaussian density 
with respect to the Lebesgue measures by generalized densities. It should be borne 
in mind that the generalized densities of Gaussian measures are defined only up 
to multiplication by a positive number, so that the above method for extending 
formulas to the infinite dimensional case applies only to formulas invariant with 
respect to the multiplication of Gaussian densities by positive numbers. The 
following propositions can be regarded as definitions of Wigner measures and 
functions similar to those given in [4], 

Proposition 3. For any density operator T on h) an ^ for ip G £2 (Q, /i), 

the following relations hold: 

(T(p)(q) = e^K 1 ^) f f e-^ qi - q) (p( qi )e- l *( B » lquqi )v[ 

Jp Jq 

(: T<p)(q ) = e*( B * lq ' q ) f [ $ T 

Jp Jq 

xe 5 ^ p,p \fi <S> v)(dq,dp). 



The notation in the first formula means that the mapping q > Wt ^ dqi + q , dpj 
is a function, while the mapping ( dqi,dp ) 1 —> Wr ^ dqi + q ; dpj is a measure. The 
function q e q - <l ) i s a generalized density of the Gaussian measure /i, and 
p 1 y P’p) i s a generalized density of the measure v. 

Let p\ be the integral kernel of a density operator T on jC 2 (Q,p) : defined by 

(Tip)(q) = e*( B ? q ’ q ) [ e^ lq ^)p(q 1 )p 1 T (q 1 q 1 )p(dq 1 ). 

JQ 
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Proposition 4. For any cp £ A)> the following relation holds 


$r (q,p) 






1 


1 


-r, qf + -r 



Let p P be the integral kernel of a density operator T on defined by 

(' T(p)(q ) = e J ( Bfl q,q ) /g zqi,qi ) p P (q, dq\). As such, p P is a function of a 

point with respect to the first argument and a measure with respect to the second 
argument 

It follows from Proposition 1 that p\(q,dqi) = f p W r ^ dgl 9 +g , dp^j . 

Setting s — r = q,s + r = q\ and using the chanmnge of variable formula, we 
obtain p P {s — r,ds + r) = J p e _1 ^ ,2r ^LFr(ds, dp), or 


Pt (q ~ 7j, dq + 0 = f p e {p,r) W T (dq, dp ); 

which means that the “measure”dp i— > Wx(dq , dp) is the inverse Fourier transform 
of the function r4 p\ (pq — dq + 0. This implies the following proposition. 

Proposition 5. Let Fe be a Hamiltonian pseudomeasure Feynman on E := 
Q x P. Then 


W T (dq, dp) = J p 2 T (q -^dq + 0 F E (dr , dp); 

Here, the integration with respect to the “measure”dg i —> Wr(dq, dp) requires 
using the so-called Kolmogorov integral @. 

2. EVOLUTION OF WIGNER FUNCTIONS AND MEASURES. 
We use the assumptions and notation of the preceding section. Suppose that, for 
each t £ R, Writ) is the Wigner measure describing the state of a quantum system 
at time t (thus in this section Wt(-) denotes a function of a real argument whose 
values are Wigner measures, while in the preceding section, the symbol Wt(-) 

5 The Kolmogorov integral is the trace on the tensor product of the space of functions on Q and the space of 
measures on Q; p\ is an element of this space (the initial definition, in which neither tensor product nor trace 
are mentioned, can be found in ED- 
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( 1 ) 


denotes a Wigner measure). Then Wt(-) satisfies the equation [5]: 

W T (t) = 2 sin 0-C u {W T (i))Y 

where a £ K, sin(aT^) is the linear operator acting on the space 7~L of E*- 
cylindrical measures on E , and conjugate to the operator sin {aC%), acting on 
the function space on E , according to 

°o 2n—1 

sin {aC%) := 'S'' —— - —^ • 

v ^ 2n-l ! n 

n =1 v 7 

(n) 

Here, C n is defined as follows: for each function T : E —> R and each n £ N, 
C^^{x) := {\Efi r), x £ E, where {T, 'H}^( a 0 := T^ n )(x)/ (g)n 'H^(x), 

TW, H (n) denote the nth derivatives of the functions in 'H. respectively, and / ,2>n 
is the n—th tensor power of the operator /, determining the symplectic structure 
on the phase space E (|5|). 

Relation (CD) implies an equation describing the evolution of the Wigner /i- 

function. To obtain it, is suffices to recall that, for any function <F : E —> R, 

the n-derivative of the product <f> n p, can be calculated by the Leibniz rule and 

that the derivatives of the Gaussian measure /i can be calculated as follows. If 

1 

h,h 1 ,h 2 ,... £ B? l Q, then [i'h = - ■) fi\ p!'h x h 2 = - {B~ x hi, h 2 ) [i + 

fi, etc. 

1 

These relations are versions of the Wick formulas. For each k £ BfrQ the 
symbol (B~ l k, ■) denotes a function defined /j- almost everywhere on Q, with the 
following properties (see Pi) 

(i) its domain is a measurable vector subspace of Q of full measure; 

(ii) this function is linear on its domain; 

1 / _1 _ 1 \ 

(hi) if x £ BfrQ, then ( B^k^x ) = ( B /c 5 13n oo J (such a function exists and 
any two functions with properties (i)-(iii) coincide /i-almost everywhere (see 

mu)). 
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For every a > 0, othe operator sin (aL^), cting on functions given on E , 
is defined by sin (aLy) (p(/i <g) v) := (sinmC^) (cpfi (8> v). Suppose also that, for 
caclit G R, is the Wigner /i-fimction describing the state of a quantum 

system at time t. 

Theorem 1. The mapping < Fr(-), taking values in the set of Wigner \i- 
functions satisfies the equation 

®T(t) = 2 sin • 

3. REDUCED EVOLUTION OF WIGNER MEASURES Let p\ n 

Py be the aforementioned integral kernels of a density operator T of a quantum 
system being the quantum version of a classical Hamiltonian system with phase 
space Ei x E 2 , where E\ = Q\ x Pi, and Eo = Q 2 x P 2 . Then, for the integral 
kernels of the reduced density operator Ti, acting on £ 2 (QpAb), i = 1,2 (here 
and in what follows, we use the natural generalizations of the above notation and 
assumptions), we have 






JQi 


the last integral is again a Kolmogorov integral. Therefore, Propositions 4 and 5 
imply the following theorem. 

Theorem 2 Let Wt and Ty be the Wigner measure and function of the 
quantum system with Hilbert space £ 2 (Q 1 x Q2, Hi ® p 2 ). Then the Wigner 
measure ITT, and the Wigner function of its subsystem with Hilbert space 
£2(^1, hi) are determined by the relations 

W Tl {dqi,dpi) = / W T (dqi,dp h dq 2 ,dp 2 ),$T(qi,Pi) 

J Q 2 XP 2 

= e K(WTi^i)+(W>i’^}) f e h({ B ^ q ^) + ( B ^ P2 )) 

J Q2XP2 

x®T(qi,Pi,q2,P2)(H2 ® zz 2 )(dg 2 , dp 2 ). 



Now, consider the models mentioned in the introduction. Throughout the rest 
of the paper, given any Hilbert space T we denote by C S {T) the set of all self- 
adjoint operators on T. 

Therefore, let 8 ^, ^i, ^ 2 , W, Wp ^2 be Hilbert spaces. We assume that 

- is the Hilbert space of a quantum control system, which we call a quantum 

plant (QF), and let W be the Hilbert space of another quantum control system, 

which we call a quantum controller (QC); suppose that hPj, and . j = 1,2 are 

the Hilbert space of parts of the QP and QC, respectively. Let M 1 := & 0 W 

be the Hilbert space of the united quantum system. Consider H^ G £' s (^), 

H<*? c ^(W), C £' S (^l 0 Wi), G We set Hf ee dback : = 

- ^— . - —-— —^ 

H^? + If ® H<^f + 0 1^20^2 + 1^!®^! 0 ^^2®^2 ^ Tfle 

G L1 s (^), G £ s (W), E^®^ g C s (& 1 (8) Wl), 2 0 W 2 ), 

are the identity operators on the corresponding spaces. The first term in the 
expression for "Hfeedback describes the evolution of an isolated QP, the second 
term describes the evolution of the isolated QC, and the last two terms describe 
the (coherent) quantum feedback. It is worth mentioning that the definition of 
^feedback is symmetric with respect to QP, QC, and the feedback. 

The more general Hamiltonian H := H^ 0 + I&> 0 H<*f + JC, where J if G 

jC s (V 0 W) ((see [ 6 ]), may describe coherent quantum control both with and 
without feedback. In particular, if JC = JPg^®^ 0 E^ 2 ®<*f 2 + ® , 

then we obtain the previous model. O 11 the other hand, if := 0 Ej 2 > 2 ®<^ 2 , 

then we obtain a model of (coherent) quantum control without feedback. 

If the QP and QC are obtained by quantizing Hamiltonian systems, then we 
can assume that, in the natural notation, Ap), = ^(Qzy O'), 

^ x Q^ 2 ,/i 1 0 /i 2 ), W = £ 2 (fe x 0 zz 2 ), j = 1,2. I 11 this 

case, the Wigner function and measure of the union of the QP and the QC are 
defined on the space x x Q^ x x Q<*? 2 , and their evolution is described 
by the equations of the second section. To obtain the dynamics of the Wigner 
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function and measure of the QP (which are defined on Q@> x x Q^> 2 ), we must 
apply Theorem 2. 

remark 1 . Obtaining the dynamics of a quantum control system (QP) 
requires finding the Hamiltonians n (c> r TO) (in appropriate classes of 
Hamiltonians). This problem is similar to the simpler problem of choosing a time 
dependent Hamilton function TO(-) on Q^ ) to which the required dynamics in 
£ 2 ( 0 ^, corresponds under the assumption = "Hi + JPi(t), where'Hi G 
J^(t) G C s (&). Although this model is not a special case of any of the 
models described above, we expect that it can be obtained as the limit of an 
appropriate sequence of these models. 

remark 2. We can extend our model, assuming that the QP interacts also 
with one more quantum system perturbing the dynamics of the control system. 
Of course, we can also assume that the source of perturbations is a part of the 

QP. 

remark 3. The approach presented in the first two sections applies directly 
to quantum systems obtained by applying Schrpdinger quantization to classical 
Hamiltonian systems. To consider more general cases, such as spin system, we 
must extend our approach by methods of superanalysis. We expect that all our 
results can be generalized to this case. 

remark 4. Feedback for classical Hamiltonian systems can be defined in a 
similar way. 

remark 5. The internal dynamics of the QP and QC in our quantum 
model with (coherent) feedback can be described in more detail. In particular, 
it can be assumed that M J = (g) 1&> 2 + (g) <S> + 1&> (g> 

1^2 “1“ ® ® + Ejq®^ ® (g> 

In the above relation, the parts of the Hamiltonian describing the QP and the 
QC and the interaction between them are again symmetric. 
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